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On  the  Hardening  Response  in  Small 


r 


Deformation  of  Metals 


by 


D,  A.  Caulk  and  P.  M.  Naghdi 


Abstract.  This  paper  is  concerned  with  a special  class  of  hardening  response 
functions  for  small  deformation  of  elastic-plastic  materials,  its  application  to 
isotropic  metals,  and  comparison  of  the  theoretical  results  with  experimental 
cyclic  stress-strain  cxrrves  for  two  different  metals.  The  theoretical  develop- 
ment is  carried  out  within  the  scope  of  the  existing  purely  mechanical  theory 
for  the  "rate-independent"  response  of  elastic-plastic  materials,  which  admits 
the  existence  of  a single  loading  function,  as  well  as  certain  accepted  idealiza- 
tions. After  summarizing  the  basic  equations  for  small  deformation,  detailed 
I attention  is  given  to  the  development  of  a special  form  of  the  hardening  response 

f , ' function,  motivated  mainly  by  the  observation  that  the  stress-strain  curves  for 

\ I uniaxial  cyclic  loading  of  a fairly  large  class  of  metals  attain  --  after  several 

j cycles  --  the  so-called  satviration  hardening.  We  exploit  this  property;  and, 

' in  the  case  of  isotropic  metals,  systematically  derive  some  restrictions  on  the 

constitutive  coefficients  in  the  loading  function  and  the  hardening 
response.  Comparison  of  the  results  with  two  sets  of  experimental  data,  obtained 
from  uniaxial  cyclic  loading  of  a 30^  stainless  steel  and  a 202k  aluminum  alloy, 
shows  good  agreement  within  the  understood  idealizations  of  the  basic  theory. 


1. 


Introduction 


By  way  of  background,  it  may  be  recalled  that  in  the  usual  development  of 
the  general  theory  of  elastic-plastic  materials,  the  rate  of  plastic  strain  and 
the  rate  of  work-hardening  are  expressed  as  linear  functions  of  either  the  rate 
of  stress  or  the  rate  of  strain.  However,  the  coefficient  response  functions  in 
these  equations,  as  well  as  the  loading  function  and  other  constitutive  response 
functions,  are  all  independent  of  the  rate  of  stress,  the  rate  of  strain  and 
time  derivatives  of  other  kinematical  ingredients.  For  this  reason,  such  a theory 
is  sometimes  called  a "rate-independent"  theory  of  plasticity,  in  contrast  to  a 
theory  which  characterizes  "rate-dependent"  behavior  by  including  rate  quantities 
in  the  loading  function  and/or  other  constitutive  response  functions.  Moreover, 
since  the  "rate-independent"  formulation  of  plasticity  theory  employs  constitu- 
tive equations  for  the  rate  of  plastic  strain  and  rate  of  work-hardening,  it  is 
also  referred  to  as  a rate-type  theory  in  contradistinction  to  a functional- 
type  theory  in  which  the  material  response  depends  on  the  entire  past  history 
of  deformation. 

The  development  of  the  present  paper  is  based  on  the  rate-type  theory  of 

elastic-plastic  materials  given  by  Green  and  Naghdi  [1,2],  which  is  constructed 

relative  to  loading  surfaces  in  stress  spacel^  We  also  utilize  the  corresponding 

fonmilation  of  the  rate-type  theory  given  by  Naghdi  and  Trapp  [3]  in  which  the 

basic  constitutive  equations  are  expressed  relative  to  loading  surfaces  in 
* 

strain  space.  Our  aim  is  to  find  an  acceptable  description  for  the  hardening 

^The  theory  of  elastic-plastic  materials  in  [1,2]  includes  thermal  effects  and 
is  developed  within  the  framework  of  a thermodynamical  theory,  while  the  develcp- 
ment  in  [3]  is  carried  out  in  the  context  of  the  mechanical  theory.  Althcjgh 
we  confine  attention  to  the  purely  mechanical  aspects  of  the  subject,  the  basic 
theory  utilized  here  can  easily  be  interpreted  in  the  context  of  the  isothermal 
theory.  The  purely  mechanical  theory  of  elastic-plastic  materials  with  large 
deformation  corresponds  to  the  general  developments  in  [1]  prior  to  thenao- 
dynamical  restrictions.  We  employ  here  a second  form  of  the  theory  of  plasticity 
relative  to  loading  surface  in  stress  space  given  in  Section  4 of  [2]. 

* 

Also,  reference  may  be  made  to  Naghdi  [4]  which  includes  a brief  summary  of 
the  main  results  in  [1,2,3],  as  well  as  a discussion  of  some  advantages  cf 
the  formulation  relative  to  loading  surfaces  in  strain  space. 
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k 


response  in  small  deformation  of  metals  within  the  scope  of  the  rate-type 
theory  of  plasticity. 

Elastic-plastic  strain  hardening  in  metals  has  received  considerable 
attention  in  the  recent  literature.  General  dissatisfaction  has  been  expressed 
with  the  more  traditional  kinematic  and  isotropic  hardening  rules,  especially 
for  describing  the  behavior  of  metals  under  cyclic  loading.  Rather  than  seek  a 
more  general  representation  of  the  hardening  response  within  the  scope  of  the 
existing  rate-type  theory  of  plasticity,  many  recent  studies  have  sought  instead 
to  introduce  modifications  in  the  basic  theory.  Before  discussing  the  nature 
of  these  modifications,  it  is  useful  to  review  briefly  certain  idealizations 
which  are  incorporated  in  the  existing  rate-type  theory  of  elastic-plastic 
materials  when  the  deformation  is  infinitesimal.  General  background  informa- 
tion concerning  these  idealizations  can  be  found  in  the  papers  of  Drucker  [5] 
and  Naghdi  [6]. 

Consider  the  uniaxial  stress-strain  curves  shown  in  Fig.  l(a)  and  l(b). 

Such  schematic  behavior  is  typical  of  ductile  metals  in  a homogeneous  deforma- 
tion. Figure  l(a)  shows  loading  in  tension  followed  by  unloading  to  zero 
stress,  followed  by  reloading  in  tension.  Figure  l(b)  is  similar,  except  that 
reloading  occurs  in  compression.  In  the  construction  of  the  basic  theory  of 
plasticity,  it  is  customary  to  idealize  the  reloading  segment  FEG  by  the  dashed 
line  FEDG.  The  loading  program  in  Fig.  1(b)  differs  from  that  in  Fig.  1(a) 
only  in  that  reloading  occurs  in  compression.  Beyond  this  fact,  the  two  schematic 
figures  are  similar,  and  the  acttial  material  response  bears  this  out:  segment  LN 
is  similar  to  EG.  It  also  seems  consistent,  therefore,  to  idealize  KLN  by  KIWN. 

In  the  same  spirit,  segments  AC  and  HJ  are  idealized  by  ABC  and  HIJ,  respectively. 

Some  of  the  idealizations  noted  above  are  incorporated  in  the  rate-type  theory 
of  plasticity  through  the  concept  of  a yield  or  a loading  surface,  which  places  a 
sharp  distinction  between  elastic  and  plastic  deformation.  If  the  yield  surface 
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is  introduced  in  the  usual  way,  then  a rate-type  theory  of  plasticity  with  I 

! 

sr^n  deformation  should  not  be  expected  to  describe  the  smooth  transition  i 

regions  such  as  EG  and  LN  in  Fig.  1(a)  and  1(b).  Recent  investigations  have 

I 

tried  to  better  represent  these  transition  regions  by  either  introducing  additional 
yield  surfaces  [7,8,9,10]  or  by  eliminating  the  yield  surface  entirely?  Eisenberg 
[12]  retains  the  idealization  EDG.(Eee  Fig.  1(a)),  but  chooses  to  describe  the 
segments  AC  and  LN  more  precisely  by  introducing  a discrete  memory  of  the  so- 
called  "last  significant  loading  event." 

At  best,  however,  the  transition  regions  described  above  are  second  order 
to  the  main  features  of  the  material  response;  and  it  seems  premature  to  intro- 
duce refinements  that  are  meant  to  improve  the  description  of  these  regions  ^ 

before  the  main  features  of  the  material  response  have  been  represented  reasonably 
well  in  the  context  of  the  basic  theory.  Certainly  the  shortcomings  of 
the  isotropic  and  kinematic  hardening  rules  are  not  entirely  due  to  idealiza- 
tions of  the  rate-type  theory.  By  the  same  token,  a detailed  account  of  the  : 

transition  regions  should  not  be  looked  upon  as  a remedy  for  poor  description  < 

i 

of  the  main  features  of  the  material  behavior.  j 

After  a brief  review  of  the  basic  constitutive  equations  for  an  elastic-  | 

i 

plastic  material  in  section  2,  we  examine  in  section  3 a linearized  version  of 
the  theory  appropriate  for  small  deformation  of  metals.  General  observations 
pertaining  to  metal  hardening,  including  saturation  hardening,  are  used  to 

( 

motivate  restrictions  and  special  assumptions  on  the  form  of  the  response  j 

^In  recent  years  some  authors  appear  to  have  placed  special  emphasis  either  cn 
the  possibility  of  developing  a theory  of  plasticity  without  introducing  a 
'yield  surface'  or  on  requiring  that  the  theory  be  capable  of  predicting  the 
existence  of  a '^ueld  surface'  rather  than  postulating  the  existence  of  a yield 
or  loading  surface  ^ initio.  As  has  been  pointed  out  already  by  Green  and 
Naghdi  [11],  such  endeavors  appear  at  present  to  be  somewhat  illusory  in  the 
following  sense:  whichever  way  the  theory  is  developed,  it  must  necessarily 
involve  some  assumptions  which  ultimately  result  in  a yield  s-jrface;  and  it 
appears  to  be  largely  a matter  of  taste  as  to  which  kind  of  assumptions  are 
preferred  at  the  outset  of  the  development  of  the  theory.  Additional  remarks 
and  some  related  references  on  this  point  can  be  found  in  [ll]. 


functions.  Isotropic  metals  are  considered  in  section  k and  the  number  of 

independent  material  coefficients  describing  the  hardening  response  is  reduced 

to  six.  These  results  include,  of  course,  the  isotropic  and  kinematic  hardening 

rules  as  special  cases.  In  section  5}  the  reduced  constitutive  equations  are 

expressed  relative  to  a loading  surface  in  strain  space  and  restrictions  on  the 

material  coefficients  are  derived.  These  restrictions  are  less  stringent  than 

those  which  can  be  deduced  via  a corresponding  development  relative  to  stress 

space.  A procedure  for  determining  the  six  material  constants  from  uniaxial 

§ ^ 

cyclic  loading  data  is  outlined  in  section  o.  In  section  7,  the  theoretical 
predictions  of  material  response  are  compared  to  some  new  and  existing  experi- 
mental data  in  uniaxial  cyclic  loading.  First,  the  cyclic  stress-strain  data 
for  30^  stainless  steel  obtained  by  Pugh  ^ [I3]  are  compared  to  predicted 

stress-strain  citrves  based  on  the  present  developments.  Next  comparison  is  made 
with  experiments  on  2024  aluminum  alloy,  which  were  conducted  at  Berkeley  by 
the  present  authors.  In  these  experiments,  cylindrical  specimens  (of  solid 
cross-section)  were  cycled  in  tension-compression  through  two  different  loading 
programs  of  two  strain  cycles  each.  For  both  the  stainless  steel  and  the 
aluminum  alloy,  comparison  of  the  theory  and  experiment  shows  good  agreement 
within  the  understood  idealizations  of  the  basic  theory. 


I 

! 

i 
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^Details  of  the  calculations  for  the  constitutive  restrictions  in  section  5 and 
for  the  determination  of  the  material  constants  in  section  6 are  provided, 
respectively,  in  Appendices  A and  B at  the  end  of  the  paper. 


4. 


2.  Notation  and  Basic  Equations. 

Let  the  motion  of  a body  be  referred  to  a fixed  system  of  rectangular 
Cartesian  coordinate  axes  and  let  a typical  particle  of  the  body  occupy  the 
position  in  some  fixed  reference  configuration.  Further,  let  x^  designate 
the  position  of  the  particle  in  the  present  configuration  at  time  t.  Then,  a 
motion  of  the  body  is  defined  by 


X.  = X. (X„,t) 

1 1 IS. 


(1) 


and  the  deformation  gradient  relative  to  the  reference  position  is  denoted  by 

3Xi 

^iK  " ‘ 


(2) 


Throughout  the  paper,  indices  take  the  values  1,2,3  and  the  usual  summation 
convention  over  repeated  indices  is  employed.  Lower  case  indices  are  associated 
with  the  spatial  coordinates  x^  and  upper  case  indices  refer  to  the  material 
coordinates  . The  Lagrangian  strain  tensor  e^g  is  defined  by 


®KL  ^iK^iL  ’ 


(3; 


where  c,„  is  the  Cauchy-Green  meas\rre  of  deformation  and  is  the  Kronecker 
delta. 

Vve  summarize  here  the  main  constitutive  equations  of  the  purely  mechanical 

theory  of  elastic-plastic  materials  contained  in  the  work  of  Green  and  Naghdi 

[1,2].  In  addition  to  the  strain  tensor  a't  each  material  point  we  admit 

the  existence  of  (i)  a plastic  strain  specified  by  a symmetric  second-order 

tensor^  e^  =e?.T(X,,,t)  which  has  the  same  invariance  properties  as  e ; (ii)  a 

scalar-valued  f’unction  K=tc(Xj^,t)  as  a measure  of  work-hardening;  and  (iii;  a 

scalar-valued  yield  (or  loading)  function  f (S|^,e^^^, k) , where  s^^  is  the 

^As  shown  by  Naghdi  and  Trapp  [l4],  symmetry  of  the  plastic  strain  follows 
as  one  result  of  a work  inequality  postulated  over  a closed  cycle  of  homogeneous 
deformation  in  strain  space. 


5. 


symmetric  Piola-Kirchhoff  stress  tensor.  For  fixed  values  of  e,„,  and  K,  the 

MN 

equation 


defines  a closed  orientable  surface  of  dimension  five  in  the  six-dimensional 
Euclidean  space  of  the  symmetric  stress  s . We  assume  the  loading  function  f 
to  ce  continuously  differentiable  with  respect  to  its  arguments. 

The  constitutive  equations  for  the  rate  of  plastic  strain  and  the  rate  of 
work-hardening,  relative  to  the  loading  function  f,  can  be  expressed  in  the  fonn 


when  f < 0 , 


= 

KL 


0 when  f = 0 and  f < 0 , 

0 when  f = 0 and  = 0 , 

A A 

X&Yj/  when  f = 0 and  f > 0 , 


where  a superposed  dot  stands  for  the  material  time  derivative  holding  X fixed 
A 

and  where  f is  defined  by 


Also,  the  scalar  \ is  restricted  to  ass'jme  positive  values  and  h,,,  aud 

KIj  ^rJj 

are  symmetric  second  order  tensors.  Each  of  the  latter  three  quantities  is  a 
function  of  the  list  of  variables 

The  four  conditions  involving  f and  f on  the  right-hand  side  of  (5)  are  called 


6. 


the  loading  criteria . Using  the  conventional  terminology,  in  the  order  listed 
they  are:  the  state  of  elastic  deformation,  unloading  from  an  elastic-plastic 


srate,  neutral  loading  and  loading  from  an  elastic-plastic  state. 

In  the  purely  mechanical  theory,  it  remains  to  introduce  the  constitutive 
equation  for  the  stress  and  this  is  specified  by 


with  the  added  restriction  that 

A A 

,,, 

' ^KL  ■ ' 

■.•■'e  r.cte  that  the  constitutive  equation  (9)^  holds  during  loading,  as  well  as 
d'uring  'unloading  and  neutral  loading,  and  that  the  condition  (lO)  ensures  that 
the  stress  is  derivable  from  a potential.  We  assume  that  for  fixed  values  of 
e:_  and  K the  expression  (9)2_  invertible  in  the  form 

A 

^KL  = 

The  preceding  summary  of  the  basic  equations  for  an  elastic-plastic 
material  began  with  equation  (4)  and  hence  is  based  on  the  concept  of  a yield 
E'urface  and  loading  criteria  in  stress  space.  Following  Na^di  and  Trapp 
it  is  possible  to  provide  an  equivalent  formulation  of  the  basic  equations 
relative  to  loading  surfaces  in  strain  space  and  this  has  certain  advantages 
which  are  discussed  in  [Si-  However,  since  a formulation  relative  to  loading 
s-urfaoes  in  stress  space  is  useful  fcr  certain  interpretations,  we  use  this 
fcrtfuiation  in  sections  3 and  4 to  motivate  special  assumptions  on  the  hardening 
response  function  and  the  loading  f-unction  f.  In  section  5j  the  results  cf 
sections  3 and  4 are  recast  in  their  corresponding  form  relative  to  loading 
S'urfaces  in  strain  space  because  in  that  form  the  derived  restrictions  cn  the 


constitutive  coefficients  are  less  stringent. 


Althoiigh  the  above  constitutive  equations  are  fully  general  in  the  context 
of  the  purely  mechanical  theory,  it  is  possible  to  conceive  of  a plausible 
physical  hypothesis  that  will  place  restrictions  on  their  form.  Starting  with 
an  assumption  for  the  work  done  by  body  forces  and  surface  tractions  in  a 
closed  cycle  of  homogeneous  deformation,  Naghdi  and  Trapp  [l4,15]  have  shown. 


among  other  results,  that  during  loading  the  plastic  strain  rate  must  satisfy 

A 


de 


aic 


MN 


A 

Ss 

Sf  MN 


Y 5 0 , 


(12) 


where  Y is  a nonnegative  scalar  function  of  the  variables  (8)  and  the  stress 
rate.  Also  recall  that  by  expressing  the  stress  response  in  the  alternative 
form  (see  [l4,l6]) 


®KL  ’ 


(13) 


the  relation  (12)  may  be  rewritten  as 


'KL  ds 


(!!a  + h . v 
™ \ p aK  MN  ^ as 


Ki  ae; 


(14) 


MN 


KL 


It  is  clear  ffom  (l4)  that  if  the  stress  response  in  (13)  depends  only  on 

that  as„-/ae^  = 0 and  as^,./aiC  = 0,  then  the  plastic  strain  rate  is 
MIN  MIN  iS-Li  MIN  I\L 

directed  along  the  normal  to  the  loading  surface  in  stress  space. 

In  the  remainder  of  this  paper  we  confine  attention  to  a discussion  of 

the  elastic-plastic  constitutive  equations  for  infinitesimal  strain.  Hence,  we 

assume  that  the  displacement, the  displacement  gradient,  the  plastic  strain  and 

the  stress  (when  expressed  in  a suitable  nondimensional  form),  as  well  as  their 

space  and  time  derivatives,  are  all  small  of  order  e.  Consistent  with  these,  the 

loading  function  f,  the  work-hardening  parameter  ic  and  its  time  derivative  tc 

must  be  at  least  of  order  € (see  the  discussion  in  section  3)-  In  terms  of  the 

components  of  the  displacement  vector  u. , the  infinitesimal  strain  tensor  e 

A KLi 

is  defined  by 


8. 


(15) 


i^iA'^A  ’ ^KL  ' -O(e)  » 


where  a comma  stands  for  partial  differentiation  with  respect  to  the  material 
coordinates  X^,  which  to  the  order  e are  coincident  with  x^.  The  infinitesimal 


elastic  strain  tensor  e is  defined  by 

K.L 


e p 

KL  KL  KL 


(16) 


We  also  recall  here  that  in  the  linearized  version  of  the  theory  the  distinc- 
tion among  the  various  representations  for  the  stress,  namely  the  nonsymmetric 
Piola-?Circhhoff , the  symmetric  Piola-Kirchhoff  and  the  Cauchy  stress,  disappears. 

In  the  next  section  we  discuss  a form  of  the  preceding  constitutive  equa- 
tions appropriate  for  metallic  materials  undergoing  small  deformations . 
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3-  Work-hardenins  of  Metals  in  Small  Deformation. 

Because  many  ideas  in  plasticity  theory  are  generalizations  of  ols-  rvations 
made  of  material  behavior  in  one  dimension,  it  is  useful  to  first  examine  the 
response  of  a typical  metal  under  a program  of  uniaxial  cyclic  loading  between 
constant  strain  limits.  Consider  the  loading  process  illustrated  by  the  stress- 
strain  curve  in  Fig.  2.  Loading  begins  in  tension  from  zero  strain  (point  O)  on 
a sample  in  the  virgin  state.  As  long  as  the  stress  remains  in  the  initial 
elastic  range,  the  stress  response  is  observed  to  be  linear  in  the  strain. 

After  the  initial  yield  point  (labeled  point  A)  is  reached,  the  plastic  strain 
increases  continuously  from  zero  until  unloading  takes  place  at  point  B.  During 
unloading  from  B,  the  stress  response  is  again  linear  in  the  strain  and  provided 
that  the  strain  at  B is  small , the  line  BC  is  usually  observed  to  be  parallel  to 
the  stress-strain'  curve  OA  in  the  initial  elastic  range. 

Motivated  by  these  observations,  we  assume  a special  form  for  the  response 
function  in  (13)  and  specify  the  constitutive  equation  for  the  stress  by 


®KL  ^KLMN^®MN‘®MN^  ^KLMN^MN  ’ 

where  the  coefficients  L are  all  constants.  From  (10)  and  the  symmetry  of 

s„,  , e.„,  and  e^„,  the  coefficients  L,,.,,,,  possess  the  symmetries 
KL  MN  RS’  KIWN 

^KLMN  " ’"LKMN  " ^KLNM  " H®KL  ' 


(17) 


(16) 


The  assumption  (1?)  also  provides  the  usual  interpretation  of  the  plastic 
strain  as  the  value  of  the  strain  at  a given  point  when  the  stress  there  is 
locally  reduced  to  zero. 

As  the  process  of  uniaxial  deformation  continues  from  point  C in  Fig.  2, 
the  material  yields  in  compression  at  point  D and  plastic  deformation  occurs 
continuously  to  point  E,  at  which  unloading  again  takes  place.  From  that 
point,  the  process  of  cyclic  loading  between  constant  strain  limits  is  repeated 


10. 


and  in  most  metals  hardening  (or  softening)  is  observed  as  the  stress  level 
during  each  succeeding  strain  cycle  increases  (or  decreases);  in  this  connection, 
a hardening  response  is  shown  in  Fig.  2.  For  most  common  metals,  the  process 
of  cyclic  loading  leads  to  a limiting  periodic  response  in  which  the  stress- 
strain  c’jrve  of  each  succeeding  cycle  is  the  same  (see  Fig.  2).  This  phenomenon 
is  sometLm.eE  called  saturation  hardening,  and  a material  is  said  to  saturate 
when  it  reaches  this  limiting  behavior.  How  rapidly  saturation  is  reached 
during  cyclic  loading  is  a property  of  the  metal. 

Based  upon  these  observations  of  the  typical  behavior  of  metals  under 
cyclic  loading,  we  assume  that  there  is  some  value  of  the  work-hardening 
parameter  v,  say  at  which  k is  zero.  Mathematically,  this  assumption  is 
expressed  as 


lim  K = lim  = 0 . (19, 

k-k:  K-K 

E S 

A sufficient  condition  for  the  restriction  (19)  to  be  satisfied  is  that  h.,,.  be 

fvi-l 

homogeneous  in  K-kg.  In  particular,  we  assume  that  h^^  is  homogeneous  of  degree 

* 

one  in  and  is  specified  by 


'^KL  -k  ’ 

o s 

where  is  the  initial  value  of  *C  and  h^^^  is  a function  of  and 
the  special  case  in  which  hj^j^  = combination  of  (6)  and  (20)  gives 


(2C 


In 


k-K 


IP 


^k  -k  ^“KL^KL 
o s 


(21; 


Thus,  in  this  special  case,  the  rate  of  work-hardening  K is  proportional  tc  the 


Since  the  experimental  observations  that  motivate  this  assiomption  do  not  depend 
on  the  deformation  being  small,  the  special  form  (20)  of  the  hardening  response 
function  may  also  be  viewed  in  the  context  of  finite  deformation. 


11. 


rate  of  nonrecoverable  work  in  the  material.  Since  a notion  of  this  kind 

is  appealing  on  physical  grounds,  it  is  desirable  to  have  the  linearized  version 

of  (6)  include  (21)  as  a special  case.  Hence,  we  assume  that  h^,  is  order  e and 

sufficiently  smooth  so  that  it  may  be  expressed  as  a linear  function  of  s and 

KL 

e^  • This  leads  us  to  write 

^KL  " \lMN^MN  ^ ^ 

where  the  constants  and  possess  obvious  symmetries. 

The  foregoing  assumptions  for  the  stress  response  (1?)  and  the  rate  of  work- 
hardening (20)  with  h^T  given  by  (22)  are  motivated  by  observations  of  the  behavior 
of  metals  in  one  dimension.  From  the  same  information,  however,  it  is  not  immediate- 
ly clear  what  may  be  an  appropriate  representation  for  the  loading  function  in  the 
infinitesimal  theory.  In  the  present  paper,  we  asstune  that  f is  sufficiently 
smooth  so  that  it  may  be  represented  as  a polynomial  in  e^  and  K.  In 

addition,  we  note  that  for  (4)  to  represent  a closed  surface  in  stress  space, 
f must  be  at  least  quadratic  in  the  stress.  It  follows  that  f must  be  at  least 
second  order  in  e.  And  whatever  the  order  of  f,  all  terms  in  its  polynomial 
representation  must  be  the  same  order  in  e.  Since  the  ratio  (k-K  )/(k  -K  ) is 

S OS 

order  1,  it  follows  from  (6),  (20)  and  (22)  that  K and  hence  also  k are  order  e . 
Thus,  the  loading  function  must  satisfy  the  relation 

f(esKLjee^,e  tc)  = e f(Sj^,e^,K)  (23) 

when 

f = 0(e'^)  , n £ 2 . (24) 

Because  for  each  value  of  the  integer  n the  loading  function  has  a different 
general  representation,  there  may  be  any  number  of  possible  forms  of  the  elastic- 
plastic  constitutive  equations,  each  corresponding  to  a different  n,  which  are 


12. 


compati'ble  with  the  assxjmption  of  small  strain.  Examples  of  two  commonly  used 
yield  functions  that  correspond  to  different  orders  of  f are  the  Tresca  and 
von  Mises  yield  conditions.  Expressed  as  a polynomial  in  s , the  Tresca  yield 
condition  corresponds  to  n = 6 while  the  von  Mises  condition  corresponds  to  n=2. 

Because  we  wish  to  describe  the  behavior  of  metals  and  since  the  von  Mises 
yield  condition  has  been  experimentally  shown  to  be  a suitable  representation 
for  the  initial  yield  surface  in  isotropic  metals,  we  assume  f to  be  second 
order  in  € . In  this  case  the  general  representation  for  f is  (see  Ref.  [l], 

Eq.  (10.11)) 


f = B 


.+  E, 


.eL  eP„  +F. 


. - »c 


'KLMN“KL“MN  “KLMN'KL’MN  * KLMN"KL''MN  ’ 


(25) 


where  the  constant  coefficients  and  have  obvious  symmetries 

KLMN  KLMN  ICLMN 

and  without  loss  in  generality  we  have  taken  the  coefficient  of  K to  be  1. 

Naghdi  and  Trapp  [15]  have  shown  that  for  a stress  response  in  the  form  (I7), 
convexity  of  the  loading  surface  in  stress  space  follows  from  s work  inequality. 
This  requirement  may  place  certain  restrictions  on  the  coefficients  in  (25). 

When  the  stress  is  specified  in  the  form  (l?))  the  relation  (l4)  for  the 
rate  of  plastic  strain  reduces  to  (see  Ref.  [15]}  Eq.  (24)) 


■P  _ 

®KL  ^ as 


, Y > 0 . 


(26) 


KL 


Since  f>0  during  loading  and  \>0  it  is  possible  to  identify 


A 

Y = If 


(27; 


and  a comparison  of  (5)4  with  (27)  yields 


B - 

Pkl  asj^  • 


(2=) 


For  example,  see  Taylor  and  Quinney  [17] j Naghdi  et  al.  [I8],  or  Bertsch  and 
Findley  [19]* 


13. 


r ’ ^ 

The  constitutive  equations  developed  in  this  section  are  appropriate  for 

i 

small  deformations  of  a general  anisotropic  metal.  The  case  of  an  initially 
isotropic  metal  is  discussed  in  the  next  section. 


4.  Isotropic  Metals. 

We  now  consider  a special  case  of  the  infinites icial  theory  appropriate  for 
metals  which  are  initially  isotropic.  Before  doing  this,  it  is  convenient  to 
introduce  some  additional  notation.  Let  s,„  and  be  decomposed  in  the  form 

lyj  ivij 

®KL  " ® ^ 3 ®MM  ’ 

(2! 

P “Po  P “P  1 P 

®KL  ® ^KL  ® " 3 ^MM  ’ 

where  and  are  the  deviatoric  parts  of  the  stress  and  the  plastic  strain, 
KL  KL 

respectively,  and  s denotes  the  mean  normal  stress . 

For  an  initially  isotropic  material,  reduces  to  an  isotropic  foiirth 

order  tensor  and  the  stress  response  (l?),  expressed  in  inverted  form,  is  given 


®KL  E ®KL  ■ E ®MM^KL  ’ 

where  E and  v are  Young’s  modulus  and  Poisson's  ratio,  respectively.  Similarly, 
the  coefficients  in  the  expression  (22)  for  h,^  reduce  to  fourth  order  isotropic 

lyj 

tensors  and  hence  this  response  function  can  be  expressed  in  the  form 

^KL  ^ “i®MM^KL‘^^2®MM^KL'^“3^KL‘^“4®L,  ’ 

where  and  are  constants.  In  addition,  the  loading  function  is 

restricted  to  be  of  the  form 

^ P2^®KL®KL^  ^4^®KL®KL^ 


where  the  coefficients  p^,...,3g  are  all  constants. 

In  the  well-known  experiments  of  Bridgeman  [20],  a large  number  of  metals 
were  subjected  to  tensile  tests  under  varying  amounts  of  hydrostatic  pressure. 
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He  observed  that  the  initial  yield  point  of  each  netal  was  not  altered  by  the  presence 


of  the  added  pressure  and  that  the  character  of  the  stress-strain  curve  d'jring  subse- 
quent plastic  deformation  remained  unchanged.  In  the  present  paper,  we  ase  the 
first  of  xhese  observations  to  place  a restriction  on  the  form  of  the  loading 
function  (32). 


Prior  to  any  plastic  deformation,  =0  and  (32)  reduces  to 

j\Li 

where  we  have  used  the  decomposition  (29)^.  Motivated  by  the  results  of 


(23) 


Sridgeman's  experiments,  we  assume  that  before  any  plastic  deformation  has 
occurred,  the  yield  function  is  independent  of  the  mean  normal  stress.  Note 
that  this  restriction  is  confined  to  the  form  of  the  initial  yield  function  and 
no  assxmiption  is  made  regarding  subsequent  plastic  deformation.  From  this 
assumption  and  (33)  we  conclude  that 


4 [f(s_^^,0,^^)]  =0  or  3Pt  + Pp  = 0 


3s 


and  hence  (32)  may  be  written  in  the  form 


1 '"2 


(3i+) 


+ a C s 

m mv 


* - * ’ 


(35) 


where  by  properly  adjusting  we  have  taken  ^^  = 1.  When  the  loading  function 
is  specified  in  the  form  (35);  the  corresponding  expression  for  the  plastic 
strain  rate  can  be  computed  from  (26)  and  is  given  by 


Contracting  this  expression  on  the  indices  K,L  and  using  {29)^,  we  obtain 

rP 


(36) 


^ = Y(3p^+Pg)e' 


(37) 


16. 


r 


"I 


If  we  assume  that  during  loading  y is  a continuous  function  of  time,  the  unique 
solution  of  (37)  which  is  compatible  with  a zero  initial  value  for  e'^  is  the 
trivial  solution 


“P  P 
6 = 

MM 


(38) 


Hence  the  plastic  volume  change  is  identically  zero.  This  result  follows 


directly  from  (3^)^  (26)  and  need  not  be  assumed  independently. 


The  expression  (35)  for  f simplifies  as  a result  of  (38),  and  the  rate  of 
work -hardening  may  now  be  expressed  as 

K-tC 


K = ( 


s xr 


)^irr 


K -K  KL’KL 
o s 


(39) 


Hence,  we  need  only  retain  the  deviatoric  part  of  h^^.  In  summary  then,  the 
loading  function  and  hardening  response  are  given  by  the  reduced  forms 


KL'KL 


'KL'^KL 


(40) 


^ ^ K -K  '■^’^KL  ’ 


o s 


where  we  have  set 


■^6  " “ ’ “3  " ^ ’ P4  = ^ ’ 


»4  = 


(41) 


.npart  from  the  term  a(Y?TY^)j  which  is  independent  of  the  stress  and  hence 

i\L  IU_i 


effectively  plays  the  same  role  as  K in  determining  the  form  of  the  yield 


s'urface  in  stress  space,  the  form  of  the  loading  function  (40)^  is  the  same  as 


one  proposed  previously  by  Edelman  and  Drucker  [22].  As  explained  in  their 


paper,  a yield  ftinction  in  the  form  (40)^  exhibits  a Bauschinger  effect  and 


reduces  to  the  von  Mises  yield  condition  in  the  absence  of  plastic  strain, 
V.'e  close  this  section  by  identifying  two  well-known  hardening  miles  , 


V.e  emphasize  that  the  restriction  (38)  is  not  an  essential  ingredient  of 


plasticity  theory  even  for  small  deformation.  It  may  be  imposed  by  sracial 
assumption  or,  as  in  the  present  development,  it  follows  from  a particular 
derivation.  In  this  connection,  see  also  [21]. 
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namely  the  isotropic  and  the  kinematic  hardening  rules,  as  special  cases  of  the 
* 

above  res’J.ts.  First,  when 


T1  = a = a = 0 and  in  the  limit  as  k , (42) 

s 

the  expressions  (40)  reduce  to 

1,2 

^ '^KL'^KL'  ^ ^ ■ (^3) 

This  is  a form  of  isotropic  hardening.  Next,  if  we  set 

T\  = ^ = 0 , ka  = , (44) 

the  expressions  (4o)  become 

f - (tj^  - 2 ' 2 ~ ) . (45) 

This  is  the  usual  form  of  kinematic  hardening. 


I • 

i 

1 ^ 

[ An  account  of  these  hardening  rules  may  be  fo\jnd  in  [6]. 
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5 . Restrictions  on  the  coefficients  in  the  yield  function  and  hardsnin/  r espouse. 


So  far  we  have  used  the  elastic-plastic  constitutive  equations  formulated 
relative  to  a loading  surface  in  stress  space  because  this  form  made  it  more 
convenient  to  appeal  to  certain  physical  observations  of  metal  behavior.  However, 
since  these  constitutive  equations  have  certain  limitations,  we  consider  further 
developments  in  the  context  of  strain  space.  Before  recording  the  corresponding 


forms  of  (40)  relative  to  a loading  surface  in  strain  space,  however,  we  need  to 
1 y‘- 


recall  fr."n  [3]  certain  results  in  terms  of  the  list  of  variables  (9)^- 


For  fixed  values  of  e^,.  and  K,  the  loading  surface  in  strain  space  is 

KL 


specified  by 


g( 


(46) 


IP 


Relative  to  g,  the  constitutive  equations  for  e^^  and  K are  expressed  as 


4P  = 

®KL 


0 


when  g < 0 , 


0 


_ A 
^KL' 


A 

when  g = 0 and  g S 0 , 


(47) 


A 

XP^g  when  g = 0 and  g > 0 , 


and 


(4S; 


where 


g = -2£-  4 
« KL 


The  positive  scalar  function  X and  the  symmetric  second  order  tensor  functions 
and  m^  depend  on 


We  do  not  elaborate  here  on  these  limitations  which  are  discussed  by  Naghdi 
and  Trapp  [ 3 ] • 
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J 


Following  [3]j  we  assume  that  g is  obtained  from  f through  the  expression 


and  that  m^j^  is  obtained  from  h^^  through  the  same  substitution  in  the 
arguments  of 

To  ensure  that  succeeding  values  of  strain  remain  on  the  yield  siorface 
, . • A ^ 

(51)  during  loading,  we  must  have  g=0  whenever  g>0.  This  condition,  together 
with  (46)  and  (49),  leads  to 

p ^ “kl  Ik ^ ^ 

A condition  analogous  to  (52)  also  arises  in  the  stress  space  formulation. 

From  (52)  and  the  sign  of  \ follows  the  restriction 


p “kl  ^ ° 

which  holds  at  all  times  during  loading . With  the  help  of  (51),  the  corresponding 

form  of  (12)  relative  to  strain  space  is 

A A 

bskl  ^ ^ as 

p dK  “mN^^MN  " ae  ■ 

To  recast  (4c)^  in  the  corresponding  forms  relative  to  strain  space,  we 

1 5“- 

first  take  the  deviatoric  part  of  the  assmed  stress  response  (30): 


^KL  = ’ 


0 0 

where  v,^  is  the  deviatoric  part  of  e and  is  the  shear  modulus  of  elasticity. 
' KL 

Combining  (5l)j  (55)  and  (40)^^  2> 
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(56) 


KL’ ^KL’ 

2 e e 
^ ■■ 

K-K 


\l 


K -K 
0 s 


(2^.SYkl  + ^Ykl)  • 


Substitution  of  (55)  into  (54)  yields 


where  the  first  of  (56)  has  been  used  to  obtain  (57)2*  Now  let 

_7\ 

Y = 2n,\g  , 


(57) 


(58) 


so  that  from  (47)  and  (58)  we  may  identify 


►'KL 


Sy 


KL 


(59) 


The  form  of  the  elastic-plastic  constitutive  equations  derived  in  section  4 
for  metals  has  now  been  completely  determined  relative  to  strain  space.  It 
remains  to  examine  the  consequenc<-E  of  the  restriction  (53)  ■ Substitution  of 
(56)  and  (59)  into  (53)  yields  the  inequality 


2p,[[4iiY^T"“Y^T.^^  )2h.y^L  (''l"2M.Q'-2a  )y^^  ] } > 0 , 


’KL  “’KL 

where  for  convenience  we  have  put 

K-K 

_s 

. K 
o s 


’’KL' 


(60) 


K-K 


3 (^^  -^  )3  J 'Tl  (j^  )'n 

^ “ OS 


(61) 


The  restriction  (60)  holds  at  all  times  during  loading . Since  it  does  not 
involve  rate  quantities,  (60)  must  then  hold  for  all  values  of  v and  K 

ivJj  KIj 

that  satisfy  g = 0.  With  the  help  of  (56)^,  g = 0 can  be  expressed  as 

^KL^L  ■ T^'^l'^L  ’ 

where  Mj^  is  the  deviatoric  second  order  tensor  defined  by 


21. 


(62) 


"kl  “ i 

For  and  K to  satisfy  g = 0,  it  is  sufficient  that  satisfy  (62). 

Thus  is  fixed  in  niagnitude,  but  otherwise  unrestricted.  Combination  of 
(63)  and  (60)  gives 

UnM^{(p+a+V)Mj^j^+  [ |a(p+a) +^-2a]Y^^}  > 0 . (64) 

which  must  hold  for  all  and  iC,  and  all  that  satisfy  (62).  It  follows 
from  the  inequality  (64)  that"^ 


P + O'  + 4g,  > 0 (65) 

and 

(i+QrH4^)^M^M^  > [i-a(9+a) +Ti- 2a]^Y^Y^  . (66) 

Since  (65)  holds  for  all  values  of  tc  and  0<(k-K  )/(k  -K  )Sl5  the  coefficients 

S OS 

a and  p musL  satisfy 

a + 4^>0  , a+0  + 4u,>O  . (67 ) 

Substitution  of  (62)  into  (66)  gives 

2 

(p+ar+4|j,)^iC  > {(p+a+4a)^(a  ' ^)  + [ 5 a(P+Q;)  + T)  - 2a]^}Y^V^  • (68) 

Fcr  a given  material,  the  inequality  (68)  must  hold  for  all  admissible  values 
cf  the  scalars  tc  and  vStvEt  • K is  bounded  by  k and  K , and  cannot  become 

-n:bo-anded  before  the  material  fractures.  Thus,  the  range  of  admissible  values 
depends  on  the  particular  material  and  a general  development 

of  the  type  considered  here  can  proceed  no  further  than  (68). 

♦ 

For  details  of  these  calculations,  see  Appendix  A at  the  end  of  the  paper. 


rcr 


K and 
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A sufficient  condition  that  (68)  be  satisfied  for  all  K and  v?..  is  that 

’ KL 

the  right-hand  side  of  (66)  vanish  identically,  i.e., 

( jap+Tl)  + ( ia^-2a)  = 0 . (69) 


j Since  this  equation  must  hold  for  all  admissible  values  of  K,  it  follows  that 

I 

! 2 
I 4a  = a and  a3  + 2T1  = 0 . (70) 

1 Substitution  of  these  restrictions  into  (56),  _ yields  the  reduced  expressions 

! 

o s 

j These  restricted  constitutive  equations  satisfy  the  inequality  (53)  identically 

i 

1 provided  a and  3 satisfy  (67),  emphasize,  however,  that  (70)  are  merely 

1 Jil  Ly2 

sufficient  conditions  that  (68)  be  satisfied;  the  only  necessary  restrictions  on 
the  constitutive  coefficients  for  a general  material  are  (67)-, 

A procedure  completely  analogous  to  the  one  employed  in  this  section  can 
' be  used  to  derive  restrictions  on  the  coefficients  a, 3, a and  T]  in  the  stress 

space  formulation  of  the  basic  equations.  Similar  results  follow  except  that 
a and  3 satisfy  the  more  stringent  conditions 


a>0  , a+3>0  (72) 

3 
i 

and  the  full  set  of  coefficients  satisfy 

> [(3+a)2(a  - ^)  + [k(P+c.)  +Ti-  2a]}Y^YKL  ^73) 

for  all  admissible  values  of  k and  The  restrictions  (72)  and  (73) 

differ  from  (67),  and  (68)  because  the  conditions  for  loading  in  each  formula- 
tion  are  not  equivalent.  To  see  this,  combine  (5l)j  (55)  and  (40)^  to  obtain 


Note  that  the  conditions  (70)^  o also  sufficient  that  (73)  be  satisfied 
for  all  K and  yBt  > jnst  as  in  ine  strain  space  formulation. 

ivL 
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6.  Determination  of  material  coefficients:  uniaxial  cyclic  loading. 
In  general,  there  are  six  material  constants 


(75) 


which  must  be  determined  by  experiment.  For  the  restricted  constitutive  equations 

(71)  , the  nvunber  of  independent  constants  is  reduced  to  four  since  in  that  case 

1 


a and  T1  are  specified  by  (70),  A procedure  for  determining  these  constants  from 
experiments  in  uniaxial  cyclic  loading  is  outlined  in  this  section. 

Consider  first  a homogeneous  extensional  deformation  in  which  the  only 
nonvanishing  component  of  the  stress  tensor  is  s^^=s(t)  say.  Then,  from  (29), 

(30),  (38)  and  the  symmetry  properties  for  isotropic  materials,  we  have 


'KL 


1 = 

0 

0 


1 

3 


0 

1 

- 3 ^ 


, / =eP  = 

’ ^KL  KL 


P 

0 


1 

' 2 ®p 
0 


0 

1 

■ 2 ®p 


(76) 


s = E(e^^-eP^)  = E(e-ep)  , (77) 

where  the  notations  e and  e^  are  introduced  for  convenience.  Corresponding  to 
the  above  homogeneous  deformation,  the  loading  functions  f and  g assume  the 
simplified  forms 

2 2 ^2 
f = 2 s -asep  + - aCp-K  , 


g = I E^(e-ep)^  - aE(e-ep)ep  +1  ae^  - K 


Also,  dtiring  loading,  the  expression  for  e^  is 


(78) 

(79) 


2 


|(p  + a + | E)e  + (Tl-|  Ea-2a)ep 


(SO) 
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where  in  obtaining  (80)  we  have  used  (52),  (57),  (58)  and  (79).  Since  the 
deformation  is  homogeneous,  we  can  obtain  the  slope  of  the  stress-strain  curve 


during  loading  by  combining  (77)  and  (80).  This  calculation  gives 


de  '■  E /r 


4 

■r  s-ae 
3 P 


(p+a)s  + 3(T)/2-a)e 


-]■ 


(81) 


The  above  expression  at  the  onset  of  initial  yield  reduces  to 


ds 

de 


. (i  + iiZi]-! 


e =0 
P 


P+a" 


(82) 


and  in  the  limit  of  saturation  hardening  the  slope  becomes 


4 ^ 

, T -r  S-ae  T 

lim  = [i  + ^ 2]-^ 

K-K  ^e  E ars-3ae 

s 


(83) 


Equations  (80)  and  (8l)  hold  for  all  metallic  materials  without 
restriction,  except  that  the  deformation  is  assumed  to  be  infinitesimal 
and  homogeneous.  If  we  impose  the  conditions  (70)  , then  these 

expressions  reduce  to 


4 • 4 

; 2.  ^ - ri  + 3 1-1 

P (F  + a + |E)  ’ de-^E  - ^ 


P + a 


(84) 


Hence,  the  slope  of  the  stress-strain  curve  derived  from  the  restricted  constitu- 


tive equations  (7l)^^2  independent  of  both  s and  e , but  is  not  constant  since 


3 depends  on  During  loading  the  slope  (84)^  decreases  from  a value  of 

4 

,-l 


ti  * 

E B + a 


(85) 


at  the  onset  of  initial  yield,  to  the  value 


[i  + 

■■E  3a-' 


(86) 


at  saturation. 

Consider  now  a uniaxial  elastic-plastic  state  specified  by 


1 


26. 


9 


s = 


^1) 


(87) 


and  suppose  that  the  material  is  unloaded  from  this  state  until  it  yields  again 
in  the  reverse  direction  at  s = s^2^-  Then,  since  both  end  ere  roots 

of  (78)  for  the  same  values  of  e^  and  k,  we  have 


..  - ^ 
(1)  (2)  " 2 


* * 
S /,  s + S 


ore 


(88) 


so  that 


* * 


a 


_ 2 !xii1!m 
3 


(89) 


p 

Experimental  data  at  a niomber  of  stress  reversals  can  be  used  in  (89)  to  obtain 
an  estimate  of  a,  and  this  value  can  then  be  used  in  (82)  to  determine  p from 
the  slope  of  the  stress-strain  curve  just  after  initial  yieldl^  Next,  a is 
determined  from  (83)  together  with  measured  slopes  of  the  stress-strain  curve 
at  various  loading  points  following  saturation.  A value  for  follows  from 
the  initial  yield  stress,  and  K is  computed  from  (78)  together  with  experimental 
values  of  stress  and  plastic  strain  at  sat\rration.  Since  T]  does  not  enter  the 
yield  function  and  it  drops  out  of  both  limiting  slopes  (82)  and  (83),  we  must 
determine  its  value  by  using  the  general  expression  (81).  Hence  we  need  to 
measure  the  slope  of  the  stress-strain  curve  at  several  loading  points  between 
initial  yield  and  saturation.  The  value  of  k at  such  a loading  point  is  computed 
from  (78).  This  value  is  then  used  in  (8I)  to  estimate  an  appropriate  value 
for  Tl. 

The  above  procedure  is  suggested  as  one  means  of  determining  the  six 


^Tien  using  this  procedure,  comments  in  the  Introduction  (Section  l)  regarding 
the  idealized  transition  from  elastic  to  plastic  deformation  should  be  kept  in 
mind.  The  gradual  transition  which  is  observed  experimentally  should  be  idealized 
by  a discontinuous  slope  in  the  stress-strain  curve  (,see  Fig.  1)  before  the  data 
is  used  to  determine  the  material  coefficients. 


27. 


■ 

material  coefficients  from  cyclic  loading  data;  by  no  means  is  the  procedure 
unique.  If  additional  experimental  information  is  available,  other  ways  of 
determining  the  coefficients  may  be  more  convenient. 


I 

I 

i 

f 


I 

I 

k 

i 
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7.  Comparison  with  experiments 

As  an  example  of  a metal  that  may  be  adequately  represented  by  the  restricted 
constitutive  equations  (7l)j  consider  the  cyclic  stress-strain  behavior  of  30^ 
stainless  steel  reported  by  Pugh  et  ad.  [13].  In  this  study,  specimens  were 
cycled  in  tension  and  compression  at  649°C  (1200°F)  between  strain  limits  of 
± 0.01  until  saturation  was  reached.  The  reported  results  [13]  are  reproduced 
here  in  Fig.  3(a).  Using  the  procedure  outlined  in  Section  6,  together  with  the 
data  of  Fig.  3(a),  one  is  able  to  choose  the  four  material  constants  in  the 


restricted  constitutive  equations  (71)-,  o 

|=4.06x10‘3  , 1=39x10 

K. 


-3 


^=0.54x10’^  , ^ = 2.46x10"^ 


(90) 


where  the  elastic  modulus  E is 


§ 


E = 123  GPa  (or  17.8x10  psi)  . 


(91) 


Theoretical  stress-strain  curves  were  determined  by  explicit  numerical  integra- 
tion of  (8o)  using  the  values  (90)  and  (91)  for  the  material  coefficients. 

Results  of  this  computation  are  shown  in  Fig.  3(t>);  and,  on  the  same  graph,  the 
corresponding  experimental  curves  from  Fig.  3(a)  are  shown  for  the  first  two 
strain  cycles.  The  comparison  for  additional  strain  cycles  is  not  exhibited 
in  Fig.  3(b)  simply  because  it  would  crowd  the  figure. 

Consistent  with  the  experimental  observations  in  [ I3] , the  theory  predicts 
that  satuiration  occurs  after  about  four  strain  cycles.  This  should  be  contrasted 
with  Eisenberg's  [ 12]  development  which  shows  at  least  eight  cycles  before 
effective  saturation.  This  occurs  in  spite  of  his  more  detailed  treatment  of 
transition  from  elastic  to  plastic  deformation. 

* 

Computations  leading  to  these  values  are  given  in  Appendix  B at  the  end  of 
the  paper . 

§9  9 / 2 

The  abbreviation  GPa  stands  for  Giga  Pascal = 10^  Pascal =10^  Newton/m  . 
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In  order  to  further  examine  the  applicability  of  the  theoretical  develop- 
ments of  the  present  paper,  we  conducted  cyclic  tension-compression  tests  on 

cylindrical  specimens  of  2024-T351  aluminum  alloy  (with  solid  cross-sections) 

* 

machined  from  plate  stock.  The  specimen  had  a gauge  length  of  2.^k  cm.  (l  in.)  and 
a diameter  of  1 cm.  (O.375  in.).  Each  specimen  was  locked  in  threaded  fixtures  and 
loaded  in  a model  TTC  SP  Instron  testing  machine  at  a nominal  crosshead  speed  of 
(C.OO5  in/min).  Strain  was  meastured  with  an  Instron  clip-on  extensometer . 

Two  different  loading  programs  were  used,  each  nonsymmetric  with  respect 
to  zero  strain.  The  results  of  the  two  experiments  are  shown  as  dashed  lines 
in  figures  4 and  5*  At  least  three  specimens  were  used  in  each  case  to  ensure 

1 

repeatability  of  the  data.  Since  relatively  high  stress  levels  were  reached 
during  the  first  loading  cycle  and  because  the  material  hardened  rather 
strongly,  the  tests  were  stopped  after  two  loading  cycles  to  avoid  buckling 
the  specimen. 

In  order  to  adequately  represent  the  observed  behavior  of  the  2024  aluminum 

alloy,  it  was  necessary  to  use  the  vmrestricted  constitutive  equations  (56),  . 

The  condition  (70)^  was  retained,  however,  for  simplicity?  Values  for  the  five 

independent  coefficients  were  selected  by  fitting  the  data  shown  by  the  dashed 

lines  in  Fig.  4.  Although  the  available  experimental  data  in  either  loading 

program  is  sufficient  to  determine  a,^  and  according  to  the  procedure  in 

section  6,  the  remaining  constants  a and  k are  difficult  to  determine  directly 

s 

since  the  experiments  did  not  proceed  to  saturation.  Consequently,  a and  tc  were 

s 

adjusted  until  the  predicted  stress-strain  curves  under  the  loading  program  in 
Fig.  4 showed  reasonable  agreement  with  the  experimental  data.  The  res\ilts  are: 

The  tests  were  carried  out  in  the  laboratories  of  the  Dept,  of  Mechanical 
Engineering  of  the  University  of  California,  Berkeley,  using  a standard 
Instron  testing  machine. 

^Recall  that  the  conditions  (TO),  „ were  only  sufficient  for  satisfaction  of 
the  inequality  (53).  ’ 
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1 = 0.07  , 1 = 0.023  , ;^=8x10'3  , 

^ = 18  X 10"^  , i ’ 

E Er 


and 


E = 69  GPa  (or  10^  psi)  . (93) 

J 

The  values  (92)  and  (93)  were  then  used  to  compute  the  stress-strain  curves  | 

corresponding  to  the  second  loading  program  shown  in  Fig.  5*  The  predicted 
stress-strain  curves  are  shown  in  both  Figs.  4 and  5 for  comparison.  It  is 
evident  that  the  nature  of  the  agreement  is  similar  in  each  case.  Also,  it  is 
interesting  to  note  that  the  reasonably  good  agreement  with  the  behavior  of 
2024  aluminum  alloy  was  achieved  at  moderate  values  for  strain,  even  though  the 
theoretical  developments  were  derived  for  small  deformation. 
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Appendix  A 

We  provide  here  details  of  the  calculations  leading  to  the  inequalities 
(65)  and  (66).  For  this  purpose,  it  is  convenient  to  regard  and  as 
vectors  in  the  five-dimensional  subspace  defined  by  the  hypersurface 
in  six-dimensional  strain  space  and  write  (62)  and  (64)  in  the  forms 


M • M = c^  , M • (aM+  by^)  > 0 , 


(Al) 


where 


= P + Q'  + 4^i  , b = '|-ci'(  p+a ) + T)  - 2 a 


KL 


(A2) 


First,  choose  M such  that  M-  y^  = 0.  Then,  (Al)  implies  that  a>0  and  the 

^ ^ <1. 

inequality  (65)  follows  immediately.  Next,  write  the  inequality  (Al)^  as 


a(M  . M)  +b(M  . ;^)  > 0 . (A3) 

The  first  term  in  (A3)  is  nonnegative,  in  view  of  (65).  Then,  by  the  Schwartz 
inequality,  namely  |M  • y^|  5 |M|ly^|,  and  the  restriction  (Al)^,  the  second  term 
in  (A3)  satisfies 

b(M.yP)  l-lb|(M.M)^(yP.yP)2  =-lbcl(yP.yP)^  , 

In  order  that  (A3)  hold  for  all  values  of  y^  and  for  all  M satisfying  (Al)  , it 

A^  X. 

D — 

is  both  necessary  and  sufficient  that  a(M.  M)  - |bc|(y^  . y^)^>0  or 

a^(M  • M)^>b^c^(y^  . y^)  . (A5) 

\ 

Next,  with  the  use  of  (Al)^,  we  obtain 
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a^(M  • M)  > b^(v^  . /) 


(A6) 


- from  which  and  (A2)  follows  the  result  (66). 
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Appendix  B 


The  material  coefficients  for  the  30^  stainless  steel  in  (90)  were  obtained 
by  the  procedure  outlined  in  section  6.  This  appendix  includes  details  of  the 
calculations • 

First,  when  (89)  is  applied  to  the  several  stress  reversals  of 
Fig.  3(a)  which  go  from  tension  to  compression,  a is  consistently  calculated 
as 

a a:  -0.5  GPa  . (Bl) 


On  the  other  hand,  the  stress  reversals  which  go  from  compression  to  tension 
consistently  give 

a 2:  1.5  GPa  . (B2) 

The  difference  between  these  two  values,  while  seemingly  large  relative  to  a, 

* * 

represents  only  small  fluctuations  in  s^^  and  s^^^  [see  Equations  (89)]  and 
may  be  due  to  a slight  material  anisotropy.  We  choose  the  average  value 

q;  = 0.5  GPa  . (B3) 


A value  for  p is  obtained  from  the  slope  of  the  stress-strain  curve  just  after 
initial  yield.  From  (82)  and  Fig.  3(a) 


ds 

de 


e =0 


P 


3.85  GPa  = [|  + 


and  hence  (33),  (b4)  and  (91)  give 


(b4) 


P = 0.477  GPa  . (B5) 

Next,  K is  computed  from  (78),  evaluated  at  e =0.  Based  on  an  initial  jneld 
o p 

stress  of  110  MPa,  we  obtain 
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tc^  = 0.0081  (GPa)^  . (B6) 


A value  fcr  iC^  follows  from  experimental  data  at  saturation.  For  example,  after 
saturation  and  just  before  unloading  from  tension,  we  have 

s = 248  MPa  , e^  = O.OO6I  . (B7) 

Substitution  of  these  variables  into  (78)  yields 

Kg  = 0.037  (GPa)^  . (b8)  i 

i 

Appropriate  normalization  of  (B3),  (B5)>  (B6)  and  (b8)  with  respect  to  the 

elastic  modulus  in  (91)  yields  the  numerical  values  given  in  (90).  | 
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Captions  for  Figitres 


Fig.  1 


Fig.  2 


Mechanical  response  of  a typical  ductile  metal  under  -uniaxial 
loading  with  theoretical  idealizations  indicated  by  dashed  lines 
(the  symbols  s and  e stand  for  one-dimensional  components  of 
stress  and  strain):  Fig.  1(a)  exhibits  loading  in  simple  tension 
followed  by  unloading  and  reloading  in  tension,  while  Fig.  1(b) 
exhibits  loading  in  simple  tension  followed  by  unloading  and 
reloading  in  compression. 


Mechanical  response  of  a typical  ductile  metal  in  cyclic  loading, 
exhibiting  saturation  hardening  and  plotted  in  the  stress-strain 


Fig.  3(a)  Experimental  stress-strain  curves  for  30^  stainless  steel  in 
cyclic  tension-compression,  ss  reported  in  Ref.  [I3]. 

Fig.  3(b)  Comparison  of  theoretical  cyclic  stress-strain  behavior  for 

304  stainless  steel  with  the  experimental  data  from  Ref.  [13]. 

The  theoretical  stress-strain  curves  ( ) are  calculated  using 

the  constitutive  coefficients  (90);  and  comparison  with  the 

experimental  data  ( ) is  shown  for  the  first  two  strain 

cycles  only,  since  curves  for  additional  cycles  would  crowd 
the  figure  [compare  with  Fig.  3(a)]. 
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Fig.  4 


Fig.  5 


Comparison  of  the  theoretically  determined  cyclic  stress -strain 

curve  ( ) for  2024- T351  aluminum  alloy  and  corresponding 

experimental  data  ( ) from  the  first  loading  program  (initially 

between  +0.02  and  -0.01  strain). 


Comparison  of  the  theoretically  determined  cyclic  stress -strain 

cvirve  ( ) for  2024-T351  aluminum  alloy  and  corresponding 

experimental  data  ( ) from  the  second  loading  program 

(initially  between  +0.01  and  -0.02  strain). 
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20.  Abstract  (continued) 

the  existence  of  a loading  function,  as  well  as  certain  accepted 

idealizations.  After  STjimnarizing  the  basic  equations  for  small  deforma- 
tion, detailed  attention  is  given  to  the  development  of  a 'Sffecial  form 
^^he  hardening  response  function,  motivated  mainly  by  the  observation 
that  the  stress-strain  curves  for  uniaxial  cyclic  loading  of  a fairly 
large  class  of  metals  attain  — after  several  cycles  --  the  so-called 
saljyation  hardening.  We  exploit  this  property;  and,  in  the  case  of 
il^ropi'c  ineiais , sysjhimatically  derive  some  restrictions  on  the 
constitutive  coefficient^,,  in  the  loading  function  and  the  hardening 
response.  Comparison  of  tls^e  results  with  two  sets  of  experimental  data, 
obtained  from  Tiniaxial  cyclic  loading  of  a 30^  stainless  steel  and  a 
2024  altuninum  alloy,  shows  gM)d  agreement  within  the  understood 
idealizations  of  the  basic  thepry. 


UNCLASSIFIED 



tKCUNITV  CLASSIFICATION  OF  THIS  FAOEr*T>*n  Dmtm  Bnfr»4) 

